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ABSTRACT 



We calculate the semi-inclusive processes B — > XsTJc and relate it to B ^ Xsip 
in next to leading order QCD, and show how it can be used to determine 
an accurate value of /^^. The exclusive modes B Krjc and B Kip also 
yield similar results, whereas B K*?7c and B K*^ can provide valuable 
information on hadronic form factors. We also estimate the branching ratios 
and the ratios of exclusive to inclusive decay modes dominated by 6 ^ srjc- 
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The b — > srjc transition offers a unique opportunity to test our under- 
standing of B meson decays. The related process b — > sV' is known to give 
the ratio for semi-inclusive decays of B, B — > ^-|- anything to exclusive decays 
B Kt/j and B K*V', in fair agreement with data [1]. Here we show that 
by taking the ratio of processes involving rjc to those involving tjj, one can 
remove the model dependence to a large extent, and have an independent 
and powerful way of determining /^^, the pseudoscalar decay constant of rjc, 
the 5*0 state of the charmonium. 

The weak Hamiltonian, corrected to the next to leading order (NLO) in 
QCD, relevant for us is given by [2,3] 



H 



AS=-1 _ 

^^=^ "V2 



.q=u,c 1=3 

where Cj are the Wilson coefficients and Oi are the operators: 



+ h.c., (1) 



01 = slY (1 - 75) bjqjr (1 - 75) Qi, 0\ = slY (1 - 75) hqjY (1 - 75) Qj, 

(2) 

O3 = s-,r (1 - 75) biJ2Vjl, (1 - 75) q'j, O4 = slY (1 - 75) bjJ27jl, (1 - 75) q'i, 

q' 9 

(3) 

05 = s-^r (1 - 75) &i e?7m (1 + 75) q'j, o, = sir (1 - 75) E?7^ (1 + 75) 

q' q' 

(4) 

The values of the Wilson coefficients at the scale /i^mb, for 

rrib = 4:.8GeV, Aj^g = 250MeV, rutop = l50GeV, (5) 

are [2]: 



ci = 1.133, C2 = -0.291, C3 - 0.015, 
C4 = -0.034, C5 = -0.010, Cq = -0.042. (6) 

We now define the matrix elements 

< O|c7^c|'0(g) >= iSn{q)9^^ < 0|c7^75c|77c(g) >= iq^fvc^ (7) 

wliere g"^ = (1.414 ± 0.083) GcV^ from ip — > e+e^ [4]. Tlic pscudoscalar 
decay constant /^^ is expected to be ~ 350 MeV, based on potential models, 
but the exact value is not known. The effective Hamiltonians in momentum 
space for the two decays are 

HeTr = ^\KsVci>\\Cr,JUq'^mh,i^ - i^MpI (9) 

where 

= C2 + C3 + C5 + (Ci + C4 + Ce) , (10) 
'c 

Cr,, = C2 + C3 - C5 + (Ci + C4 - Ce) . (11) 

We shall treat and C^^ as a phenomenological parameters, thus ab- 
sorbing in theirs definition any higher-order correction or deviation from 
factorization that may arise [1]. The value of |C^| determined from BR{B 

Xs4>) = (1 ± 0.25) X 10~^ [5], where we have corrected for the fact that 
not all if] 's are produced directly but arise from cascades of if]' 's, etc., is 

|Cv,| = 0.220 ±0.026. (12) 
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This value of \C^\ favors =2. Wc note that the ratio \Cr^^/C^\ is one in 
leading order. In the NLO, with the QCD coefficients from Ref's [2] and/or 
[3], for the 2 < A^e < 3 and = 174 ± 25 GeV the ratio is stable within 
the following values: 



I C^7C^ I = 0.89 ±0.02. (13) 

The ratio of semi-inclusive production to 77c production can now be 
easily calculated. We find that 



where 



T(h ^ sr]c 



.9^ . 



A 



1/2 



X 



r(6 ^ s^p) 

[(mg -m^^^(mg + mg)] 

[ml{ml + m^) - m^(2m^ - m^,) + mj 

- 4.0(GeV-2) 



A: 



be 



2 2 \ ^ 2 2 

ml ml) mi 



(14) 



Measurements of this hadron-model independent ratio offers a very accurate 
determination of fj^^. In ratio (14) , again only the direct 7]^ production should 
be included. Since the cascade of charmonia into r]c is extremely small, no 
correction is necessary. 

Next we consider the B — > and B — > Kric modes. Using the general 
Lorentz decomposition of the matrix element 

< K{k)\-s^,h\B{j)) >= (p + A;)^/S(g') + gMA"i(g'), (15) 



4 



we have 



1/2 



X 



1 - 



m, 



Vc fi~). 



m, 



Since = m^^, we have set /^(?7^^J//xb("^^V') ~ ^' '^^^ second term 
in Eq.l6 is [mlJm\] = -0.06, and it is neghgible 

compared with the first term in (16), in accord with the pole approximation 
as discussed in [6]. Because of that, essentially hadron-model independent, 
ratio (16) is: 

^ 2 

V{B Kric)/T{B KiP) = 14.2(GeV-2) 



(16) 



Finaly, we calculate the exclusive rates for B K*il) and B — ; 
The general Lorentz decomposition of the (B— > K*) matrix element: 



(17) 



< K*{k)\sj^{l-^,)b\B{p) >= -te^,^,{p + kf{p-k)\"{k)V 

+e{k){ml - ml.)Ai - {q ■ e{k)){p + k)f,A2 (18) 
+{q ■ e{k)){mB + mK*){q^,/q'^) [2mK*Ao - (m^ - mK*){Ai - A2)] , 

and the corresponding definitions of relevant form-factors 



(niB + rriK*) 



_ Aq (q-) / 2N _ 

^0 — 7 : r, ^o[q j — j-^ — 

[mB + mK*) ^ ^ (1 - 



y(o) 



(l-?7m?_) 



^(0) 



q^/ml-) 



(19) 



(20) 
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(niB - rriK*) 



^i(O) 

- q^/ml+) 



(21) 



A, 



I \ V' ^2 j — Y^ 

[mB + mK*) ^ ' (1- 



^2(0) 



(1 - q^lm\^) 



(22) 



^o(O) 



^(0) 



mB - rriK* 



2mK* 2mK* 
gives the hadron-model dependent ratio 



^2(0), 



T{B K*ip) 



/ niB + rriK* \ 



94' 



3/2 



l^lol 



X 



+ 



2\vr + 



+ 



m. 



The ratio (24) can be represented as 



m 



K* 



nil 



1 - 



m 



K* 



(23) 



(24) 



m. 



1 - 



r(S ^ K*r]c)/T{B K*iP) = R{GeV-^) 



Jr]c /-< 



(25) 



The factor R depends critically on the form factors. We evaluate it in three 
models for illustration. The first is the relativistic quark model, Model 1, 
[7]. For Model 2, we use a parametrization based on effective chiral the- 
ory for mesons with flavor and spin symmetries of Heavy Quark Effective 
Theory [HQET] [8,9]. This method contains errors that follow from semilep- 
tonic D-meson decays used as inputs [10]. The Model 3 is a variation [11], 
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where A2{0) is set equal to zero and V(0) is renormalized from the value in 
Ref. [7] . Last model gives a much better fit to experimental measurements of 
^l/^t+l ratio, where I^LiT) refers to the longitudinal (transvers) polarization 
of mesons m. B ^ K*'^! decay. This ratio is also very sensitive to the form 
factors. All models use the single pole fits for V and A^, (i=0,l,2) form- 
factors. The masses of excited meson states, corresponding to the sh current, 
are: mo- = 5.38, mi- = 5.43, mo+ = 5.89 and mi+ = 5.82, GeV's. In Table 
1 we summarise the parameters (V,Aj) and the factor R for S — > K*{i/j,ric) 
decays, while the ratios Tl/I't+l, ^t/^t+l are given for B — > K*il) decay. 



Table 1 

The (V,Ai), R, Vl/Vt+l and Tt/Vt^l in Models 1, 2, 3 and 1', 2'. 



Model 
[Ref] 


V(0) 


^(0) 


^i(O) 


^2(0) 


R 

{GeV-^) 




Tt/Tt+l 


1[7] 


0.369 


0.321 


0.328 


0.331 


1.01 


0.573 


0.427 


2 [8,9] 


0.61 


0.20 


0.20 


0.20 


0.67 


0.366 


0.634 


3 [11] 


0.656 


1.133 


0.328 


0.0 


5.49 


0.728 


0.272 


1' 


0.37 


1.13 


0.33 


0.0 


6.0 


0.80 


0.20 


2' 


0.61 


0.69 


0.20 


0.0 


4.7 


0.63 


0.37 
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The experimental value of Tl/Tt^l is: 

VL{B^K*i^) J 0-80 ±0-051^8 ,[12] 
Vt+l{B ^ K*ij) I Q gg ^ O.lOiH^ .[13] 



(26) 



Note that if the condition ^2(0) = is put by hand in Models 1 and 2, we 
find Tl/Tt+l — 0.80 and 0.63, respectively. In Table 1 we call these models 
1' and 2'. Measurement of D ^ K*lh' decay [10] favors ^2(0) = 0, and 
HQET relates {B K*) to {D K*) transition [11]. 

As a by-product of the above analysis and using recent CLEO II experi- 
mental results on the B — > K*'4> and B — > K-^ decays [12]: 

BR{B- K*-^p) = (0.178 ± 0.051 ± 0.023) x 10"^ (27) 

Si?(S° ^ X*V) = (0.169 ± 0.031 ± 0.018) X 10"^ (28) 

BR{B- K-^p) = (0.110 ± 0.015 ± 0.009) x 10"^ (29) 

BR{B^ XV) = (0.075 ± 0.024 ± 0.008) x 10-^ (30) 
from the expressions, 



T{B^X,7],) A.{)V{B^K*'n, 



T{B Xs^p) R T{B -> K*tlj) ' 

T{B^ XsVc) ^ 4.0 TjB ^ Kri,) 
r{B Xsij) 14.2 r{B Kijj) ' 

we estimate the ratios of exclusive to inclusive decay modes: 



(31) 
(32) 



r{B- K*-r]c)/T{B X,r]c) = R{QM5 ± 0.018), (33) 

r(£° ^ V) /r(5 ^ XsVc) = R{0M2 ± 0.014), (34) 



T{B- K-r]c)/T{B X,r]c) = 0.391 ± 0.116, (35) 
r(S° ^ K%c) /T{B XsTic) = 0.266 ± 0.112. (36) 

The above estimates reflect the experimental errors from (27) to (30). 
TheoreticaUy, the ratios in (33) and (34) have to be equal, as also (35) and 
(36). If one uses the values of R from models based on A2(0)=0, then r{B — > 
K*ric) is = 25% of the inclusive rate. This seems to be surprisingly large 
because it would imply that B — > Krjc and B — > K*r]c predominantely 
saturate transition B Xgrjc ! 

To estimate branching ratios for B — > K*rjc, B — > K 77c and B — > Xgrjc 
decays one has to know the pseudoscalar decay constant /^^. Theoretically, 
like the value of g^, the value of /^^ can be related to the wave function of 
the ^-state of the charmonium at the origin: 

4 = 12m^|V^(0)|^ (37) 




Without QCD corrections the expressions (37,38) gives /^^ = 350 MeV. The 
QCD corrections are significant but approximately cancel in the ratio [14], 
and one expects /^^ to be close to the above value. 

A non-perturbative estimate of /^^ based on QCD sum rules [15] could 
be more reliable. Saturating vacuum polarization of pseudoscalar density 
evaluated to 0{as) with m^^, m^/ and using n=4,5 moments, which minimize 
the continuum contributions, we find 



= (300 ± 50)MeV. 
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(39) 



The major source of error above is not as{m,^J = 0.212 ± 0.021, but the 
current quark mass of charm which we take as rric — (1.25 ± 0.03)GeV^. 

Using the central value of /^^ = 300 MeV, and taking the ratio |C^^/C^| 
= 0.89, we estimate branching ratios for B — > K*r]c, K r]c and Xgrjc decays: 



BR{B- K*-r)c) = i?(0.127 ± 0.040) x 10'^, (40) 

BR{B'^ K*^r]c) = i?(0.121 ± 0.026) x 10~^ (41) 

BR{B- K-7]c) = (1.113 ± 0.177) x 10-^ (42) 

Si?(S° ^ K%^) = (0.759 ± 0.256) x 10-^ (43) 

BR{B X.rie) = (2.851 ± 0.713) x 10"^ (44) 



To conclude, we have shown a very accurate technique of measuring /j,^. 
The consistency of this value obtained from semi-inclusive decays and exclu- 
sive decays will check factorization. Once /^^ is known from the measure- 
ments of B — > XsTjc, this can be used to see if the form-factor ratio (24) 
is consistent with models of those decays. The measurement of B — > K*rjc 
probes the spin-0 part of the axial form factor, and provides a useful check 
on the model building. Finaly, we estimated branching ratios and the ratios 
of exclusive to inclusive decays for the most interesting modes. 

This work was supported in part by the U.S. Department of Energy under 
Contract DE-FG06-85ER40224, and by the Croatian Ministry of Research 
under Contract 1-03-199. 
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